2011 Fermi Symposium, Roma., May. 9-12 



1 



Three-dimensional Two-Layer Outer Gap Model: the Third Peak of Vela Pulsar 
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Department of Physics, University of Hony Kony, Pokfulam Road, Hony Kony 

We extend the two-dimensional two-layer outer gap model to a three-dimensional geometry and 
use it to study the high-energy emission of the Vela pulsar. We apply this three-dimensional two- 
layer model to the Vela pulsar and compare the model light curves, the phase-averaged spectrum 
and the phase-resolved spectra with the recent Fermi observations, which also reveals the existence 
of the third peak between two main peaks. The phase position of the third peak moves with the 
photon energy, which cannot be explained by the geometry of magnetic field structure and the 
caustic effect of the photon propagation. We suggest that the existence of the third peak and its 
energy dependent movement results from the azimuthal structure of the outer gap. 



I. INTRODUCTION 

The FermVs observation of the Vela Pulsar [l[ 
shows that, in its light curve, as the energy increases, 
a third peak appears at the trailing part of the first 
peak and shifts towards the second peak, which makes 
this known 7-ray pulsar before Fermi's era more spe- 
cial. 

We want to use our two-layer Outer Gap model 0] 
to explain the strange third peak and its movements. 
Our two-dimensional two-layer model has explained 
the phase averaged spectra of the mature pulsars, in- 
cluding the Vela pulsar, in the first catalogue of 7-ray 
pulsars of Fermi [3j. However, to study the pulse pro- 
file, a three-dimensional model is needed. Therefore, 
we extend the two-layer model to a three-dimensional 
magnetic field. Here we present our simulations of 
the energy dependent light curves and the explana- 
tion of the moving third peak of Vela Pulsar, by three- 
dimensional two-layer Outer Gap model. 



II. THE THREE-DIMENSIONAL 
TWO-LAYER OUTER GAP 

In this study, the three-dimensional rotating vac- 
uum dipole field is adopted. The 7-ray emissions come 
from the region with strong accelerating electric field, 
which extends above the last-open filed lines and be- 
tween the null charge surface of the Goldreich- Julian 
charge density [4j and the light cylinder. This region 
has a two-layer structure. As shown in Figure [1] in the 
trans-field direction of the magnetic field, the outer 
gap can be divided into two parts: 

1 the main acceleration region at the lower part 
of the outer gap, where the charge density is 
assumed to be ~ 10 % of the Goldreich- Julian 
value and a strong electric field is accelerating 
the particles to emit GeV photons via the cur- 
vature process, 

2 the screening region around the upper bound- 
ary, where the growth of the main acceleration 



region in the trans-field direction is stopped by 
the pair-creation processes. 
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FIG. 1: The simplified distribution of the charge density 
(solid line) and the corresponding accelerating potential 
(dot-dashed line) of the two-layer outer gap. 

And we use a simple step function to approximate 
the distribution of the charge density in the trans- 
field direction in the poloidal plane (the plane where 
the field lines have same polar angle p )as 

, , s _ / pi(x,cj)p), if < z < hi(x,(t> p ) 
p{X ' Z ^ p) ~\p 2 (x,(l> p ), if fcifo^,) < z < h 2 {x,4> p ) 

(1) 

where x, z and 4> p represent coordinates along the 
magnetic field line, the height measured from the last- 
open field line, and the azimuthal direction. In addi- 
tion, hi and h 2 represent the thickness of the main 
accretion region and the total gap thickness, respec- 
tively. 

Instead of solving the real three-dimensional Pois- 
son equation, we divide the gap into many equal divi- 
sions in the azimuthal direction and regard each slice 
of the gap as a two-dimensional two-layer gap. For 
each division, the solution of the potential in our two- 
dimensional study is applied. The shape of the 7-ray 
spectrum of the division is determined by three pa- 
rameters: /, fractional gap thickness, h\/h 2 , ratio of 
the thicknesses of the primary and whole region, and 
pi, the number density in the main acceleration re- 
gion. 



eConf C110509 



2 



201 1 Fermi Symposium, Roma., May. 9-12 



We define the gap fraction / measured on the stellar 
surface as @, 



_ w 

J — • 



(2) 



where i? s is the stellar radius, and r p (4> p ) is the po- 
lar cap radius. The accelerating electric field £|| is 
proportional to f 2 . 

The accelerating electric is caused by the devia- 
tion of the charge density from the Goldreich- Julian 
charge density, so we introduce a parameter g to rep- 
resent this deviation, the g and p satisfy the relation, 
p(x,z,<p p )-pGj(x,(f>p) ~ g(z,(f> p )pGj(x,(f>p)- For each 
slice of the gap, the deviations of the charge density 
in the two region, — g\ and (72, and h\jhn satisfy 



III. 



(^) 2 = 1+^. 
hi 92 



SIMULATION OF THE ENERGY 
DEPENDENT LIGHT CURVES 



(3) 



The first step is to determine the viewing angle 
and the inclination angle by making the 'geometry 
determined' light curve, whose peaks are due to the 
caustic effect and almost do not dependent on the en- 
ergy. Our two-dimensional study shows that the frac- 
tional size of the gap is small (around 0.16), there- 
fore, in this study, we trace the magnetic field lines of 
a = 1 — > 0.93. The factor a is used to represent the 
magnetic field lines at a given layer, where a — 1 and 
a = represent the last-open field lines and magnetic 
axis, respectively. 

The emission direction is calculated in observer's 
frame 5]. The curvature photon is assumed to be 
emitted in the direction of the particle motion, which 
can be described as an combination of the motion 
along the magnetic field line and the drift motion, 



v p B/B + rxfl, 



(4) 



where v p is calculated from the condition that \v\ = c. 
The polar angle to the rotation axis £ of the emission 
direction and the pulse phase ip are calculated from 



cos C, = v z /v 
ip = — cos^ivx/vxy) ~ r ■ v/Ri 



(5) 



If a point in a field line satisfies |£ — /3\ < ip(r), the 
radiation of that point can be seen by the observer, 
where /? is the viewing angle and (p(r) is the solid angle 
of the radiation. 

We find that the set of (a, /3, a rn i„) = 
(57°, 80°, 0.935) reproduces right peak separation of 
the two peaks. Figure [5] shows the skymap of the emit- 
ted photons from the magnetic surface of a = 0.95 
with a = 57° and (3 = 80°. 
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FIG. 2: The skymap of the radiations outside the null 
charge surface with inclination angle a = 57°. The z-axis 
is the pulse phase and the y-axis is the direction of the 
radiation. The dashed line is the viewing angle, which is 
chosen as f3 — 80° . 



Then we calculate the 7-ray spectrum to find proper 
values of the three parameters, 1 — gi, h\jhi and /. 
We find that the values of these three parameters used 
in the two-dimensional study can provide a good spec- 
trum. Here we assume the distance to the pulsar is 
325 pc. 

At last, we calculate the phase resolved spectra un- 
der the parameters found above, and integrate the 
phase resolved spectra to obtain the energy dependent 
light curves. The number of the photons measured at 
pulse phases between tpi and ip2 is calculated from 



N^(Ei,E 2 ,ipi,ip2) oc 



E-2 



E 1 



Ftot(E,ipi <tp< tp2)dE, 
(6) 

where the F to t is the phase resolved spectrum. 

The calculated pulse profile under constant values 
of the three parameters in azimuthal direction can not 
explain the existence of the observed third peak. The 
fact that the simple caustic model can not explain 
the phase shift of the observed third peak, forces us 
to consider more complex structure of the emission 
region. 

By the definition of / in equation we may 
choose the form of the azimuthal distribution of / as, 



f(<f> P ) = 



C 



(7) 



where the C 



0.18r™ Qa: , 



is the maximum value 



of the polar cap radius, and the factor 0.18 is chosen 
by fitting the phase averaged-spectrum. 

It is expected that as the null charge surface is closer 
to the stellar surface, the number density of the X- 
ray photons increases in the gap and the screening 
region becomes thinner due to higher pair creation 
rate. Therefore, we assume the formula of the ratio as 



hi 



P ) - B 1 + B 2 



l/r n uii((f>p) - 



max 
null 



l/, 



,min 
null 



l/r; 



max 
null 



(8) 
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TABLE I: The effects of the distributions of the three parameters 



fiib) 


hi/h 2 {<t> p ) 


Pl{<t>p) 


The Bump 

O.lGeV < E < lGeV 


The Third Peak 

IGeV < E < 8GeV 


The Third Peak 

8GeV <E< 20GeV 


0.2 


0.927 


0.05 


no 


no 


no 


0.2 


Equation [5] 


0.05 


phase ~ 0.2 - 0.3 


no 


no 


0.2 


0.927 


Equation 0& [9 


no 


phase ~ 0.26 


phase ~ 0.26 


Equation [7 


0.927 


0.05 


no 


no 


phase ~ 0.4 


Equation [7 


Equation [8] 


Equation 0& 


phase ~ 0.2 - 0.3 


phase ~ 0.26 


phase ~ 0.33 



By fitting the spectral shape, we obtain B\ = 0.89 
and B-2 = 0.09. 

With the two-dimensional two-layer model, we 
found that the phase-averaged spectra for most of 
the 7-ray pulsars can be reproduced by the averaged 
charge density of p = [hipi + (h 2 - hi)p 2 ]/h 2 ~ 
0.5. In three-dimensional magnetosphere, the actual 
Ej_(4> P ) can make the particles in different <^ p -cell to 
drift into other t^-cell via the effect of E± x B. As a 
result, the actual average density is 



P(<P P ) = Po 



f{4> p + &4> p ) 
f(4> P ) ' 



(9) 



where the p(<fr p ) = N \4>p) / f '{4>p) IS use d and p — 0.5 
is the averaged charge density without drift motion. 
Using the relationship between g\, g 2 and h\jh% given 
by equation ([3]) , the distribution of 1 — g\ can be ob- 
tained. 

We expected that the particles' displacement due 
to the drift motion becomes more important on the 
magnetic field line that has a smaller radial distance 
to the null charge surface, because the particles run 
longer distance in the outer gap. To take into this 
effect, we assume the formula of displacement as 



= F 



l/r nu ii((f>p) - 1/r. 



max 
null 



l/r. 



null 



i/r: 



max 
null 



(10) 



where F is a fitting parameter, which is chosen as -28°. 

Table Q] summaries the effects of the distributions of 
the parameters on the energy dependent light curves. 
Figure [3] shows the pulse profiles calculated by taking 
into account the azimuthal distributions /, h\/h% and 
1 — gx . The corresponding phase averaged spectrum 
is shown in Figure 2] Figure [S] is the intensity map 
in the pulse phase and the energy plane. The color 
represents the scaled number of the photons at the 
certain interval of the pulse phase. We find a moving 
third peak from the results. 

We also fit the calculated phase resolved spectra 
with power law plus exponential cut-off form. Fig- 
ure [6] shows the cut-off energy and the photon index, 
as functions of the pulse phase. 
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FIG. 3: The energy dependent light curves (histogram) 
with the distributions of 1 — g\, /u/fe and /, provided by 
the three-dimensional outer gap. The solid lines are the 
observed light curves from Fermi-LAT Q. 



IV. THE REASON FOR THE THIRD PEAK 

The reason for the third peak and its shift is as 
follows. In the light curves, when energy lower than 
1 GeV, the distributions of the thickness ratio h\jh 2 
makes a third-peak structure at ~ 0.2 pulse phase, 
because h\jh 2 affects the emissivity in the screening 
region, which mainly produces the curvature photons 
of the energy less than 1 GeV. In the energy bands 
higher than 1 GeV, on the other hand, the azimuthal 
distributions of the fractional thickness, /, and the 
number density, 1 — g±, are more important, and they 
produce the third peak at ~ 0.3 — 0.35 pulse phase by 
providing higher electric field and more particles, re- 
spectively. Consequently, the differences in the stand- 
ing phases of the third peak due to the distributions 
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FIG. 4: The phase averaged spectrum with the distribu- 
tions of 1 — gi, hxjhi and /, comparing with the observed 
data (circle) from Fermi-LAT [1[ 




FIG. 6: The cutoff energies and photon indices of the 
phase-resolved spectra of different pulse phases, comparing 
with the observed data (circle) 

and j^{4>p), when pi and hi/h 2 get close to and 
1, respectively. This makes the effects of azimuthal 
structure of the gap obvious. 4) The inclination angle 
also determines the azimuthal structure of the gap. 
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FIG. 5: Intensity map in the pulse phase and energy plane. 
The color represents the percentage of the number of the 
photons of certain interval of pulse phase in the total num- 
ber of photons of certain interval of energy. 



V. CONCLUSION 

The third peak of Vela Pulsar can not be explained 
by the caustic effect determined by the magnetic field 
structure. This peak is due to the azimuthal distri- 
butions of the fractional gap thickness (/), the ratio 
of the thicknesses of the primary and whole region 
(/11//12), and the number density in the main acceler- 
ation region (pi ) . The distributions of p\ and / make 
third-peak-like structure in the bridge region of light 
curve above 1 GeV, while the distribution of h\jhn 
makes a bump in the bridge region of the light curves 
below 1 GeV. The phases of the third peaks caused by 
the azimuthal distributions of h\jhi, p\ and / are dif- 
ferent from each other. Consequently, the differences 
in the phases produce the shift of the combined third 
peak with the photon energy. 



of hi/h 2 , f and pi(orl — gi) produce the shift of the 
third peak with the photon energy. 

Why is the Vela Pulsar so special? This is because: 
1) It has a thin gap to make sure that its observable ra- 
diation contains both the contributions of the two lay- 
ers. 2) Its 7-ray radiation is dominated by curvature 
photons that keep the effects of azimuthal structure of 
the gap visible. 3) Two-dimensional study shows that 
the shape of the spectrum is very sensitive to pi(4> p ) 



Acknowledgments 

We thank L. Zhang for useful discussion. This 
work is supported by a GRF grant of the the Hong 
Kong SAR Government entitled "Gamma-ray Pul- 
sars" (HKU700911P). 



[1] Abdo A. A. et al, 2010, ApJ, 713, 154 
[2] Wang Y., Takata J. & Cheng K. S., 2010, ApJ, 720, 
178 

[3] Abdo A. A. et al., 2010, ApJS, 187, 460 

[4] Goldreich P. & Julian W.H., 1969, ApJ, 157, 869 



[5] Takata J., Chang H.-K. & Cheng K.S., 2007, ApJ, 656, 
1044 

[6] Yadigaroglu I. A., 1997, Ph.D. thesis, Standford Uinv. 



eConf C110509 



